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We present a new model, based on colour rearrangements, which at the same time 
can describe both hidden and open charm production in S-meson decays. The model 
is successfully compared to both inclusive decays, such as B — > J/ipX and B — > D S X, 
as well as exclusive ones, such as B — > J/ipK^*' and B — > D^D^K. It also gives 
CN ! a good description of the momentum distribution of direct J/^'s, especially in the 

J> | low-momentum region, which earlier has been claimed as a possible signal for new 

' exotic states. 
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INTRODUCTION 



A proper understanding of the confinement phenomenon in quantum chromodynamincs 
(QCD), describing the transition from the perturbatively calculable parton level to the ex- 
perimentally observable hadron level, is still missing. In order to get a better understanding 
of the hadronisation process one is therefore lead to constructing models, such as the Lund 
string fragmentation model and then compare these to data. Of special interest in 
these types of models is the treatment of the colour quantum number and the associated 
colour flux. Typically, the planar approximation is used for this, which is valid in the large 
Nc — > oo limit, leading to a good description of inclusive event properties at high energy 
colliders. 

At the same time, there is also a large class of so called hard diffractive processes in both 
ep and pp collisions which cannot be described by the planar approximation. The signifying 
feature of this type of processes is that the final state particles are divided into two or 
more colour singlet systems, separated by large rapidity gaps, and that at least one of these 
^ ■ systems has the properties of a hard partonic interaction, such as jets. These events are in 
accordance with the predictions of the model introduced [2j based on pomeron exchange. 
Although further developments of this model is successful in describing rapidity gap data, 
it implies different descriptions of diffractive and non- diffractive events without a smooth 
transition in between. 

In order to order to remedy this situation and get a model that can describe both inclusive 
and hard diffractive processes, the Soft Colour Interaction (SCI) model was introduced [§]. In 
short, this model is based on the assumption that the colour flux from the hard perturbative 
interaction is modified through interactions with the colour background field represented by 
the remnants of the incoming hadrons. In the simplest version of the model there is only one 
new additional parameter describing this interaction, namely the probability for such a colour 
exchange. These colour rearrangements were added to the Lund Monte Carlo programs 
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Lepto [4] for deep inelastic scattering (DIS) and Pythia [5| for hadron-hadron collisions, 
where they give rise to events having regions in phase space where no string is stretched, 
and therefore no hadrons are being produced. These rapidity gap events are classified as 
diffractive and this simple model essentially reproduces all data on diffractive hard scattering 
in both ep [3| and pp (sf collisions. Of course, events without such rapidity gaps are also 
produced corresponding to ordinary inclusive events. The SCI colour rearrangements can 
also turn a colour octet cc or bb pair into a singlet, leading to production of charmonium 
and bottomonium states in basic agreement with data from pp j3] and pA and ttA jsf . 

This phenomenological success of the SCI model ^ indicates that it captures some very 
essential QCD dynamics. It is therefore interesting that recent developments on QCD rescat- 
tering theory 1(| provides a basis for this model 11]. Rescatterings of a hard-scattered 



parton on the spectator system cannot be gauged away and do contribute at leading twist. 
In DIS such rescatterings of the struck quark via 1,2. . . gluons are summed in the Wilson line 
used in the definition of the parton density functions, which thereby absorb these rescatter- 
ings effects when fitted to inclusive DIS data. However, for less inclusive observables that 
depend on the colour structure in the event, these rescatterings are important and the SCI 
model is a phenomenological model to account for their effect. 

There are also several extensions of the model. The difference in the potential energy of 
various string configurations can be included 13] and the momentum transfer in the colour 
exchange can be modelled The model has also been successfully extended to describe 
jet quenching in a quark gluon plasma [13 . 

The rational behind the SCI model is to learn more about the non-perturbative dynamics 
by starting from a well defined perturbative state. One such example, which is the subject 
of this paper, is the production of hidden and open charm in hadronic -B-decays. Thanks to 
the 5-factories there is now a wealth of detailed data from the BaBar and Belle experiments, 
which can be used as a testing ground for the ideas behind the SCI model. In a I?-meson 
decay the hard scale is given by the 6-quark mass, m& ~ 5 GeV, and the decay products 
then interact with the remaining soft part of the 5-meson. Following earlier applications of 
the SCI model, we are aiming at formulating a model which can describe both hidden and 
open charm production. 

This paper is organized as follows. We start in section 2 by reviewing earlier models for 
hidden and open charm production in B-decays, with emphasis on charmonium production, 
and how they compare with data. In section 3 we then present our model on the parton level 
and the transition to the observable hadron level. The resulting model is then compared to 
existing data in section 4 and finally section 5 contains the conclusions. 



II. EARLIER MODELS 

Charmonium production in .B-decays has a long history. In a naive version of the so called 
colour singlet model (CSM) [l4| one notes that there is a probability 1/9 that the cc pair in 
a b — > ccs decay is in a colour singlet state and combines this with the J/ ip wave-function at 
the origin giving a rate which is in reasonable agreement with data. However, in a proper 
treatment one also has to take into account that the Hamiltonian describing the 6-decay 



should be the effective one, where the H^-boson (and t-quark) has been integrated out [15 
This gives an additional factor ~ 15 suppression of the colour singlet state compared to the 
colour octet one To make things even worse, a strict NLO calculation of the colour 



singlet rate becomes negative 17] unless one also includes 0{afy corrections to the colour 
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octet channel and thus modifies the perturbative expansion. Even so, the rate obtained 
is about a factor ten below data. In addition the CSM cannot be used to calculate the 
production of P-wave charmonium states (Xc)- 

A theoretically more sound description of charmonium production in S-decays is provided 
by the so called colour octet model (COM), based on non-relativistic QCD HI ], in which 



the decay is factorised into two parts. First, a non-relativistic colour singlet or octet cc- 
pair is produced in a given spectroscopic state, 2S+1 Lj, where S, L and J are the spin, 
orbital angular momentum and total angular momentum respectively, and then this state 
is transformed into a charmonium hadron through the possible emission of soft gluons in 
order to get a colour singlet system. The latter process is described by non-perturbative 
matrix-elements which in principle can be fitted from data. 

If one only takes into account the colour singlet and octet 3 Si states and uses data from 
the J/if> decay width and the production of J /if) at large transverse momenta at the Tevatron 
to fit the non-perturbative matrix elements, then the COM predicts a rate for direct J /if) 
production which is about half of the observed value 19| and similarly for if)' . However, in a 
proper non-relativistic expansion the colour octet 1 5 , and 3 Pj, J = 0, 1, 2 states also need to 
be included [l6| . Unfortunately, the corresponding non-perturbative matrix-elements cannot 
be fitted independently so it is not possible to get a prediction for these rates. Once the 
J I if) rate has been fitted it can be used [l6| to get a prediction for the inclusive t] c rate, for 
which there is yet no data. Finally one can fit the branching ratio B — > Xc2X and get a 
prediction for B — > XciX which turns out to be about a factor two below the observation. 

Given that the J / if) rate is fitted to data, the COM can then be used to calculate the J/ if) 
momentum distribution by assuming that the 6-quark decays through a two body process 
b — > (cc)q, taking into account smearing from the boost to the B (or T(AS)) rest system [2(| 
as well as a non-perturbative "shape-function" which resums the emissions of multiple soft 
gluons Overall this gives a good description of the data on direct J /if) production when 
combined with a model for the J/if>K(*) contributions except for small momenta where it is 



much below the data (22j. Several suggestions have been made for how to describe this low- 



momentum region, including enhanced baryon pair production [23j , hybrid mesons [2J, |25 
and hidden charm diquark-antidiquark bound states [26l |. 

As already alluded to our model aims to describe both open and hidden charm production. 
Whereas the production of charmonium states in 5-decays have been studied extensively, 
there has been less attention devoted to B — > DD S and B — > DDK decays. There are some 
models, which all are more or less based on the factorization hypothesis and to varying degree 
make use of heav y qu ark symmetry. For example, for B — > DD S decays we have the simpler 
pole models (27l . 28 ] which later have been refined using heavy quark symmetry. 



Another example is given by [31(, which is an extension of the so called Isgur-Scora-Grinstein- 
Wise model [32| to non-leptonic decays, also assuming factorisation. Typically all of these 
models are able to describe the branching ratios of the specific decay processes they are 
studying (see for example (3l[) but they cannot be used to predict hidden charm production. 
Finally, there also exists models for B — > DDK decays assuming that they proceed via an 



intermediate DD sJ state [33j, |34|, which however are not so successful at describing data. 



III. OUR MODEL 



The model has three main ingredients: The internal 5-meson dynamics with the b- 
quark decay giving a partonic final state, the soft colour interactions which modifies the 
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colour structure of the event and finally hadronisation using the Lund string model amended 
with special treatment of systems of small invariant mass and the mapping onto discrete 
charmonium states. It is particularly important to have a properly devised and tuned 
hadronisation model in order to describe the exclusive few-body final states in .B-meson 
decays that is investigated here. In this section we define and discuss these details of our 
model. 



A. .B-meson dynamics and 6-quark decay 

The .B-meson can be viewed as a non-relativistic 6-quark surrounded by a non- 
perturbative hadronic system represented by a spectator quark to account for the light 
valence quark as well as seaquarks and gluon contributions. This is the basis for the AC- 
CMM model j35|, which describes the internal dynamics of the .B-meson. In the .B-meson 
rest system the three-momentum p b of the 6-quark is spherically symmetric and given by 
the normalised Gaussian 

having a width given by the parameter pp. The other model parameter is the mass m sp of 
the spectator. The decaying 6-quark is not a final state parton and its mass m b is allowed to 
vary dynamically as given by energy-momentum conservation, M\ = (p b +Psp) 2 , resulting 
in 

m 2 b = M*+m%-2M By Jpl + m%. (2) 

The weak decay of the 6-quark is illustrated in Fig. [TJ which defines the four-momenta 
in the process, where in the b rest frame Q = (m^, 0). The differential decay rate in the b 
rest frame is then 



d r v „ = k Wc1 ^ SHQ - Pq - Ph - fh )JL± E l^ 2J g; )3 4^,3 2E %, ■ 

(3) 

where ifh/si are ^-factors for the hadronic and semi-leptonic decays, respectively, that will 
be fitted to data to get the correct normalization. The leading order spin averaged squared 
matrix element for the 6-quark decay can be expressed as [36] 



\Y,\M\ 2 = ^G F \V hh \ 2 \V qb \\p b -p h ){ Pq -p h , 



2^—' — — ' - — ^ {k2 _ M 2 w)2 + T 2 wM 2 w (4) 

In our Monte Carlo model we thus start by generating the three-momentum pf, of the 
6-quark from Eq. (TO , which is then used in Eq. (j2J) to get the dynamical 6-quark mass m b . 
The momenta of the decay products from the 6-quark decay are then generated according 
to the differential decay rate in Eqs. (13141) in the 6-quark rest frame and boosted to the 
.B-meson rest frame. 



B. Colour structure 

From the calculation of the two colour configurations in the b — ► ccs decay using the 
effective theory [ill , we know that the cc colour singlet fraction is suppressed with about 



FIG. 1: Feynman diagram of the weak 6-quark decay into a lighter quark, q, and a fermion pair, 
fi and f2- Q is the momentum of the incoming 6-quark, k is the momentum transferred by the W 
and p the momenta of final state particles. 
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FIG. 2: Feynman diagram showing the colour string connection with normal colour connection (I) 
and with reconnected colour (II). The double dashed lines indicate the colour string connections. 

a factor 100 compared to the colour octet one. Therefore, we will in the following simply 
assume that all the parton level decay products from the 5-meson are in the colour config- 
uration represented by diagram I in Fig. [21 where cq is one colour singlet system and cs the 
other as indicated in the figure. 

In order to model the colour suppressed mode where the cc-pair forms a colour singlet, 
as illustrated in diagram II in Fig. [2J we use the Soft Colour Interaction (SCI) model 0]. 
This model is based on the assumption that partons emerging from some hard, perturbative 
process interact softly with the colour background field provided by the spectators of the 
initial hadron as they propagate through it on a Fermi length scale. In these soft processes, 
the small momentum transfers are not important and can be neglected, at least in a first ap- 
proximation. Instead, it is the colour exchange that is important, since it changes the colour 
string topology of the event and thereby affects the hadronisation giving a different hadronic 
final state. The model uses an explicit mechanism where colour-anticolour, corresponding 
to a non-perturbative gluon, can be exchanged between partons and remnants. This should 
be a natural part of the process when bare partons are dressed into non-perturbative ones 
and the confining colour flux tubes (strings) is formed between them. 

The SCI model has been added to the Lund Monte Carlo programs Lepto [4] for deep 
inelastic scattering (DIS) and Pythia j^] and for hadron-hadron collisions. The hard parton 
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level interactions are given by standard perturbative matrix elements and parton showers, 
which are not altered by the softer non-perturbative effects occurring on a longer space- 
time scale. The probability for a SCI in terms of the exchange a soft gluon within any 
pair of a parton and a spectator remnant cannot be calculated and is therefore taken as 
a constant given by a phenomenological parameter P, which is the only parameter of the 
model. As mentioned in the Introduction, this model has support from QCD rescattering 
theory and is very successful in describing data on hard diffraction, i.e. rapidity gap events, 
and charmonium production. 

Applying the SCI model on I?-meson decays means that the perturbative partons from 
the 6-quark decay will undergo soft non-perturbative interactions with the background colour 
field in the 5-meson represented by the spectator quark. Again, the parameter P specifies 
the probability for such a colour-anticolour exchange between any of the partons with the 
spectator. Since there are here only 4 partons, including the spectator, there are only two 
possible string configurations 1 as shown in Fig. [2J Starting from configuration I, a soft gluon 
exchange switches the system to configuration II. A second such colour exchange switches 
the system back to I again. Thus, increasing the colour exchange probability P too much 
will not favour configuration H. This switch-back effect is extreme in this 4-parton state. 
A similar effect appears for the rapidity gap rate in high energy ep and pp collisions, where 
additional colour exchanges may switch back from a gap topology to a no-gap topology, but 
is less pronounced due to the presence of more partons. As shown below, the fit to S-meson 
decays gives the value P ~ 0.15, whereas the rapidity gap rate is not strongly dependent on 
P and the model is stable for P ~ 0.2 — 0.5 [f|. 

C. Hadronisation of small mass string-systems 

In order to be able to compare our model with data we also need to describe the transition 
from the parton level to hadrons. As starting point we use the Lund string model [l[ as 
implemented in the Pythia Monte Carlo |5j together with its special treatment of string 
systems with small invariant masses, which produce only one or two hadrons. For 5-meson 
decays the latter part of the model is most important and, as we will see below, we have 
had to introduce a more careful treatment of these small-mass systems in order to describe 
data on exclusive decay modes. In addition, to calculate the probabilities for different 
charmonium states we use the model as will be discussed in the next subsection. 

In the standard Lund string fragmentation picture hadrons are produced iteratively by 
considering pair production of quark-antiquark pairs 2 in the colour field of the string leading 
to the production of two colour singlet systems: one of which is a quark-antiquark pair which 
becomes a meson and the other which is a rest string. This process is then repeated until 
the rest string has such a low invariant mass that the procedure is terminated by producing 
either two mesons or a single one as will be discussed below. 

Normally, for a given quark anti-quark pair that is going to form a meson, the model in 
Pythia only produces, i.e. maps the pair onto, the two lowest order mesons, i.e. those with 



1 Due to the small phase space for extra emissions the parton shower does not give any additional gluons 
and strings for a large majority of decays. 

2 In order to be brief we only describe here the production of mesons. The Lund string fragmentation model 
can also describe the production of baryons, for example by production of diquark-antidiquark pairs. 
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L = 0, S = 0, J = (such as if-mesons) and L = 0, S = 1, J = 1 (such as AT*(892)) 
mesons. However, to be able to fit both inclusive, where all type of mesons contribute, and 
exclusive branching ratios, where mostly the lowest order mesons contribute, we also need 
to activate the production of Axial Vector Mesons (AVM) with L = 1, S = 0, J = 1 (such as 
i^i(1270)). The production of these mesons are controlled by the PYTHlA-parameter Pavm 
(PARJ(14)), which is the probability that a S = meson is in a L = 1 state. These mesons 
have substantially higher mass than the ordinary ones, which necessitates refinements to the 
model for their production in small mass systems. 

Consider a low mass system, with invariant mass m sys . To hadronise this system Pythia 
first tries n try (MSTJ(17)) times to make two hadrons with total mass below m sys . The two 
hadrons are given a relative momentum to conserve the invariant mass of the system. If 
the program fails to make two hadrons, one hadron is made instead and put on-shell by 
exchanging an effective gluon with some other part of the event. We have changed this 
second step to also try n try times to make a single hadron with mass below m sys and if this 
fails accept the last tried hadron even though m h > m sys . 

In the one hadron case, there are several ways to put the hadron on-shell and obey energy- 
momentum conservation. In default Pythia it is done in different ways depending on the 
situation: If the system in question is the only one left to hadronise, it exchanges momentum 
with the, already produced, final state particle which is furthest away in momentum space. 
When another unhadronised system exists two different procedures are used. If rrihad < msys, 
the four-momentum vector, Phad, is scaled down and the excess momentum is put as a gluon 
in the other system. If rrihad > m sys, the four-momentum needed is taken from the other 
system. 

We have tried a number of different detailed treatments for momentum exchange to 
unhadronised systems. It turns out that this part of the model has little influence on 
the branching ratios, but can affect the J/^-momentum distribution. The best result is 
obtained when requiring that the effective gluon has p 2 = 0. In practice this is done by 
giving the gluon a fraction x of the system three momentum and (1 — x) to the hadron, i.e. 
p g = (x\p sys \,xp sys ) and p had = (E sys - x\p sys \, (l-x)p sys ), with the condition p\ ad = m\ ad . 
Fig. O shows the mass, energy and momentum distributions of the exchanged gluon for the 
different parton system cases. It can be noted that p 2 is not always zero, because if there are 
no unhadronised systems to exchange momentum with then momentum is instead exchanged 
with another final state particle as described above. 

D. Hadronisation to charmonium states 

A special case of hadronisation of small mass systems is the mapping of colour singlet cc 
systems produced at the parton level with a continuous mass spectrum onto the discrete mass 
spectrum of charmonium states. To calculate the probability to obtain different charmonium 
states we use the model js[ , which is based on the assumption that it is more likely that a cc 
pair of given mass m cg is mapped to a charmonium state which is close in mass rather than 
to one which is further away. This can be motivated by the fact that the charmonium mass 
spectrum covers a mass range of almost 1 GeV, which is substantially larger than the few 
hundred MeV energy-momentum transfers of the soft colour interactions that may affect the 
invariant mass of the system. For example, a cc with mass just above the threshold 2m c , 
using m c = 1.35 GeV, should have a larger probability to produce a J /if) than a if)', and a 
cc close to the open charm threshold should contribute more to if)' than to J /if). 
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FIG. 3: Distribution in mass, energy and three-momentum for the effective gluon, with p 2 = 0, 
which is transferred from the specified qq-p&ir, when it hadronises into only one particle. Indicated 
in the figure are the color connections for the pair, see Fig. [2j 



Thus, the model assumes that the smearing of the cc mass due to soft interactions is 
described by the Gaussian 



(m r r — m} 2 



G sme (m c5 , m) = exp f ^ J , (5) 

V ° sme / 

where the width is a sme = 0.4 GeV. The probability that a cc pair of mass m cE forms a 
charmonium state i of mass rrii is then given by 

„/ f G sme (m C c,m)Fi(mi,m)dm , s 

l^jj G sme {m c ^m)F j {m j ,m)dm 

where i^(mj, m) is the distribution in invariant mass for a given charmonium state i. For 
our purposes it is enough to use the approximation Fi{ra^ m) = Si5(m — rrii), i.e. neglecting 
the very narrow width of charmonium states but including the relative weights Sj = 2 Jj + 1 
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FIG. 4: Probabilities used for mapping cc pairs of mass m c g onto different charmonium states and 
open charm states (DDX). 

coming from non-relativistic spin statistics. 3 The expression we use is then 

7-,/ \ srne yTfi C c, Tfii) 

Pi{m C c) = = - r. (7) 

In our study the following six charmonium states have been included, r) c , J/ip, Xoc Xic 
X2c and ip(2S). Figure H] shows the resulting probabilities for the different charmonium 
states. It includes the decrease in total probability for charmonium production due to open 
charm production and the effect from trying more than once to make a single hadron giving 
an additional suppression of heavy mesons at low m c5 . 



IV. RESULTS 

Before comparing our model with data we want to emphasize the importance of the model 
used for hadronisation. To illustrate this Fig. [5] shows the mass, energy and momentum 
distributions for the parton string-systems obtained from the matrix-element for each of the 
possible colour configurations compared to the resulting final state hadrons. As can be seen 
from the figure this mapping is far from being smooth, especially when there is only one 
hadron produced from the initial qq pair. It should also be noted that energy and momentum 
conservation implies that the energies of the sc and cq systems, and of sq and cc, add up to 
the -B-mass, as well as that their momentum distributions are pairwise the same. 



3 We do not need to include any additional suppression factor 1/n for states with higher main quantum 
number n, which was included in the earlier studies 0, Isf, since the model used then did not include 
the additional suppression of heavy mesons from trying more than once to make a single hadron with 
miiad < m sys . 
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FIG. 5: Distributions of mass, energy and momentum for partons and final state hadrons. Dashed 
lines show the distributions of the initial partons. Solid lines show the distributions of the final 
state particles originating from the parton pair, q is either a d or a u quark depending on if the 
meson is a B° or a B~ . Indicated in the figure are the color connections for the pair, see Fig. [21 



A. Normalization 

As specified in Eq. (j3J) we have included two different if-factors, one for semi-leptonic 
and one for hadronic decays, in order to get a correct normalisation and to take into account 
the difference between semi-leptonic and hadronic decays. We also be note that, since the b- 
quark mass is varying in the underlying 5-meson model and the decay width is proportional 
to the 6-mass to the fifth power, the K-factors are sensitive to the parameters of the B- 
meson model and not necessarily larger than 1. The two A"-factors have been obtained by 
simultaneously fitting the branching ratio for semileptonic decays b — ► e~X to data using 
B°/B ± -> l+X = 10.24% 0, and also fitting the total width from all the decays included 
in Tabledto the measured lifetime. The resulting A"-factors, for pp = 0.57 GeV, m sp = 0.15 



GeV and assuming a total width of 4.20 x 10"" GeV [37 



are 



K sl = 0.75, K h = 1.4 (8) 

and the resulting b branching ratios are shown given in Table [H Of special interest here is 
the value for the Br(b —>■ c s c) (marked with boldface) which corresponds to the partonic 
mode that has been used to fit the model and which is in good agreement with data, 
Br(b^ esc) = (22 ±4)% Q. 
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decay channel 


CKM- 
factors 


Simulated decay 
width (1(T 15 GeV) 


Branching ratio 
without iT-factors 


Branching ratio 
with two IC-factors 


b- 


> c e v e 






16.08% 


10.14% 


b- 


> u e~ v e 




0.75 


0.21% 


0.13% 


b-. 


• c pT 


Kb 


55.71 


15.99% 


10.08% 


b^ 


■ U\L~~ 




0.72 


0.21% 


0.13% 


b- 


> cr v T 


K b 


14.11 


4.05% 


2.55% 


b^ 


■ U T~ V T 


Kb 


0.27 


0.08% 


0.05% 


b - 


->■ c s c 




63.40 


18.20% 


22.09% 


b- 


-* c s u 


V c bV us 


6.86 


1.97% 


2.39% 


b- 


-» c d c 


VcbV cd 


3.54 


1.02% 


1.23% 


b- 


c d u 


VcbVud 


143.85 


41.29% 


50.11% 


b- 


-> u s c 


VubVcs 


1.08 


0.31% 


0.38% 


b- 


-> « s u 


VubVus 


0.09 


0.03% 


0.03% 


b- 


u d c 


VubVcd 


0.06 


0.02% 


0.02% 


b- 


-* u d u 


VubVud 


1.92 


0.55% 


0.67% 


All channels 




348.35 


100.% 


100.% 


Total fraction of c 






118.% 


122.% 



TABLE I: Decay widths and branching ratios for 6-quark decays used to fix the lf-factors. The 
semileptonic decays, marked with italic, and the total width are used to fix the factors K s \ and 
Kfr. The channel b — ► c s c is used to normalize our results. 

B. Summary of the model 

The complete new model contains 4 parameters: 

• the width pp of the Gaussian 6-quark momentum distribution in the 5-meson, 

• the mass m sp of the spectator quark in the 5-meson 

• the probability P for soft colour exchange, 

• the probability -Pavm for producing a meson with L — 1, S — 0, J — 1 

The parameter values given in Table [III are determined by fitting to S-meson decay data. 
The values of pp and m sp of the ACCMM-model are sensitive to the momentum spectrum of 
the produced J/ip, whereas P and Pavm are fitted to the branching ratios. This is done by 
simulating a large number of B decays where b — *> c s c and then calculating the branching 
ratios using the normalization described above. A x 2 is then calculated as 

2 _ (BR exp — BRdata) / Q \ 

Two sets of branching ratios have been used for the fits, one with only inclusive branching 
ratios and one with both inclusive and exclusive ones. When changing pp and m sp within 
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PdVCl TYl py 




LLlX^l Llol V c J-J1V 


Fitted to 

all 


/ / In TY1 MTTl OTltll TTl 

U 1 iy IllwiIlcIlL U.II1 


PF 


0.30 - 0.60 






0.57 GeV 


Wisp 


0.05 - 0.15 






0.15 GeV 


p 


0-1 


0.16 


0.16 




Pavm 


0-1 


0.71 


0.72 





TABLE II: The parameters of the model: the width pp of the Gaussian 6-quark momentum 
distribution in the i?-meson, the mass m sp of the spectator quark in the i?-meson the probability 
P for soft colour exchange, the probability -Pavm for producing a meson with L = 1, S = 0, J = 1. 
Their values are obtained by fitting to: the inclusive branching ratios in Fig. [6l all branching ratios, 
i.e. also including the exclusive ones in Fig. and the J ftp momentum spectrum in Fig. 



the ranges given in Table [TTl it is always possible to find a reasonable fit to the branching 
ratios. We have tried a number of different combinations of parameter values, including 



the ones used in |20|, [2JJ and the one that gives the best fit to the direct J/ijj momentum 



distribution is then used. 

As discussed in Section [Til C\ decays to mesons with L = 1 (e.g. i^i(1270)) requires to 
account for production of the much heavier axial vector mesons and we obtain Pavm ~ 0.7 
for the probability that an S = meson is in an L = 1 state. We note that this value 
is close to 3/4 as obtained by simple counting of available angular momentum states for 
L = and 1. (In default Pythia Pavm — since observables based on 'stable' hadrons are 
not sensitive to whether these mesons have been produced as intermediate states or not.) 
The inclusion of such heavier mesons requires a retuning of the Pythia parameter n try 
controlling the number of tries allowed in particle formation from small mass systems. By 
increasing from the default value n try = 2 one accounts better for the available phase space, 
and we have chosen n try = 6 as preferred when fitting to all branching ratios. Fitting only 
inclusive branching ratios would prefer a slightly higher value, nt ry = 9, but this degrades 
the fit to all branching ratios more than using nt ry = 6 does for the inclusive fit and, in 
addition, gives a worse fit to the J/ip momentum distribution. 



C. Comparison to data 

1. Inclusive branching ratios 

Fig. [6] shows the results of the fit to the inclusive branching ratios as well as the fit to 
all branching ratios compared to data. As is clear from the figure, in both cases the model 
describes the inclusive branching ratios very well. In fact, when fitting to these inclusive 
channels only, we obtain 

y 2 5 

= -f = 0-7 . (10) 
d.o.j. 7 

We also note that the difference between the two fits is very small, it is only barely visible 
for the B — ► J/i/jX and B —>■ D S X channels. 
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FIG. 6: Comparison between data [37( (red dots with error bars) of specified inclusive decay 
channels for i?-mesons and the model fitted to these data only (green down-pointing triangles) and 
to all branching ratios (blue up-pointing triangles). The lower part is for B°, the middle part for 
B~ and the top part for B° and B~ combined. Note that the B — > J/ipX channel is multiplied 
with a factor 10 to make it more visible. 

2. Exclusive branching ratios 

Fig. [7] demonstrates that also the exclusive branching ratios are quite well described 
by the model and that the difference between the two fits is again very small. In both 
cases the model essentially describes all the different channels involving D-mesons, including 
the relative strength of D and D* channels. It also gives an overall good description of 
the states that can be produced from the two colour configurations displayed in Fig. [2] 
and discussed in Section IIII Bl respectively. This is true both for the two-body decays 
shown in Fig. E](c,d) as well as for the three-body decays shown in Fig. [7](a,b) where for 
some channels both colour configuration I and H contribute as indicated in the figure. As 
a consequence the model also gives sum-rules of the type, Br(B° — > K°) ~ 

Br(B° -> D^ + D^°K-) + Br(B° -> D^°D^°K ), which, at least within errors, are in 
agreement with data. Evidently it is more demanding to describe all these exclusive final 
states which are sensitive to the non-perturbative dynamics in hadronic few-body systems 
with relatively small kinetic energy available. In view of this, it is remarkable that this 
model with only a few parameters reproduce the data so well. We also note that the rate for 
J/ipKi (1270), which is controlled by the -Pavm parameter, comes out essentially right even 
if this decay mode is not included in the fit. 

Not surprisingly, however, there are some channels where the description is not so good. 
This is mainly for channels with one of the heavier charmonium states and either a K 
or a K*, e.g. B° — > \icK*° . In our approach the two mesons are formed more or less 
independently and we therefore get about the same ratio for K° and K*° irrespectively of 
whether it is produced together with a heavier or lighter charmonium state. However, the 
data indicates that the production of mesons with small relative momenta, i.e. closer to 
threshold, such as XicK* , should be suppressed. 
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FIG. 7: Comparison between data (red dots with error bars from 38( in (a,b) and from 37j] in 
(c,d)) of specified exclusive decay channels for B° (a,c) and B~ (b,d) and the model fitted to only 
the inclusive decay channels in Fig. [6] (green down-pointing triangles) and also including these 
exclusive channels (blue up-pointing triangles). The lower parts are associated with the original 
colour string configuration I, whereas the upper parts (with upper BR scale) are for configuration 
H (cf. Fig. [2]). In (a,b) there are middle parts where both I and II contribute. 



When fitting to all channels we get x 2 /d-o.f. = 353/46 = 7.7, which may be considered 
too large for a good fit. Most of the x 2 i s > however, coming from some particular channels. 
For example, removing the B ip(2S)K channels, which gives the dominant contribution, 
results in x 2 /d.o.f. = 181/44 = 4.1. We also note that since our model is not based on 
first principles, it is not very meaningful to perform x 2 -tests of it. Instead, it is meant to 
investigate whether the SCI model, which is phenomenologically very successful in describing 
other kinds of data related to colour string-field topologies, is of relevance also in decays of 
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Direct J/v|/ momentum (GeV/c) 



FIG. 8: Momentum distribution of directly produced J /tp in B — » J/ijjX decays in the T(4S) rest 
frame. Data from 39|] and 22| compared to different models: our soft colour interaction model 
fitted to only inclusive decay branching ratios and to all branc hing ratios, and the NRQCD curve 
which is a combination of the COM [2lf] and J/tpK(*) done in 22]. 



5-mesons, and this seems indeed to be the case. 



3. The J ftp momentum distribution 

Finally, in Fig. [S] the results of our model is compared to the momentum distribution of 
directly produced J/^'s in B decays measured in the T(4S) rest frame. First of all it is clear 
from the figure that our model gives a good overall description of the data especially in the 
low-momentum region, whereas it is slightly too high in the peak region. At the same time 
it should be noted that the model has not been fitted to the normalisation of these data, 
only to the inclusive and inclusive plus exclusive branching ratios respectively. Since the 
peak in Fig. [S] is dominated by two-body decays, the model's excess here of ~ 20% indicates 
that it produces somewhat too many direct J/ip's, in particular in the J/ipK channel as is 
also indicated by figure Fig. [7(c). 

For comparison Fig. [H]also shows the results of a NRQCD based model from j22| discussed 



in Section 2. This model is a combination of the COM results [21] together with a model 
for the exclusive J/ipK and J/ipK* decays. In contrast to our model, the NRQCD based 
model is not at all able to describe the low-momentum region, which has given rise to various 
alternative explanations as already mentioned. It is also important to recognize that the 
parameter that mostly affects the overall shape of the momentum distribution, namely the 
width, Pi?, of the Gaussian momentum distribution of the 6-quark in the I?-meson, does not 
affect the low-momentum region of the spectrum. Instead this parameter is responsible for 
the smearing in the peak region. Similarly, the mass m sp of the spectator quark also does 
not affect the low-momentum region. In addition, for both of these parameters we have used 



the same values as in [20j. The main difference of our model compared to the NRQCD-based 
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ones therefore lies in the dynamical treatment of the decay products from the 6-quark. Such 
soft dynamics is particularly important to get a good description of the low-momentum 
region in Fig. [HJ and our detailed model does indeed provide an improvement. 

V. CONCLUSIONS 

In this paper we have presented a new model, which can describe both open and hidden 
charm production in 5-meson decays. The model is based on the ACCMM model for the 
internal 5-meson dynamics, the SCI model for colour rearrangements in the partonic final 
state, a model for mapping of colour singlet cc pair on to charmonium states and a new 
procedure for hadronising colour singlet systems with small invariant mass within the Lund 
string fragmentation framework. 

Using more or less standard values for the parameters of the ACCMM and SCI models: 
the width of the Gaussian momentum distribution of the 6-quark in the I?-meson, pp = 0.57 
GeV, the mass of the spectator quark, m sp = 0.15 GeV, and, the probability for a colour 
rearrangement, P = 0.16, we find overall good agreement with the data both on open 
and hidden charm production. In order for this to be possible we have had to activate 
the production of axial vector mesons in the hadronisation with a probability -Pavm = 0.7 
consistent with the number of available angular momentum states. We have also improved 
the probing of the available phase space, which is particularly constrained for processes 
close to mass thresholds, by increasing the number of times that the program tries to make 
a single hadron out of a given small mass partonic system. Related to this, we have modified 
the way that energy and momentum is exchanged when the invariant mass of such a parton 
system has to be changed in order to give the proper hadron mass. 

Our model gives a very good description of inclusive observables such as B — > J/ipX 
and B — > D S X, which shows that the basic ideas ingredients of the model are correct. In 
particular it shows that the idea of soft colour interactions also can be successfully applied 
in 5-meson decays giving a unified description both of open and hidden charm production. 
When it comes to exclusive decay modes the overall description is still good but there are 
some channels, which are not well described. The latter is especially true for B — > ip(2S)K^ 
and B — > XicK* indicating that there is something lacking in the model for mapping colour 
singlet cc pair onto charmonium states which also may be connected to the fact that these 
decay modes are the ones closest to threshold that we have considered. Last but not least, 
the model also describes the momentum distribution of direct J/ip's, including the low mass 
region where earlier models, most notably the COM, fails. This shows that there is no 
need to invoke new hadronic states such as hybrids or bound diquark anti-diquark states to 
explain this region. Instead this comes out naturally from our model as a consequence of 
the non-perturbative dynamics involved in the hadronisation process. 

The overall conclusion of our study is that it is possible to describe and understand a 
wealth of data on 5-meson decays with our relatively simple model based on the framework 
of Soft Colour Interactions, which modify the colour structure of an event and thereby the 
string-field topology leading to different hadronic final states. The fact that this SCI model 
has earlier been successful in describing a wide range of phenomena, such as rapidity gap 
events and charmonium production in both hadron-hadron and electron-proton collisions, 
and now 5-meson decays shows that it captures essential generic features of non-perturbative 
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QCD interactions. 
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